*3 * V, 4 




SEMIANNUAL REPORT FOR 


Grant No. NAG 1-500 




SEMICONDUCTOR SUPERLATTICE PHOTODETECTORS 


January 1, 1986 to June 30, 1986 


Submitted to 


Dr. Ivan Clark 

National Aeronautics and Space Administration 
Langley Research Center, Mail Stop 473 
Hampton, VA 23665 

PASA-CE-1793 50) SEMICONDUCTOR SUP ERL ATT ICE 
PHOTO DEI EC TOR S Semiannual Repcit, 1 Jan - 
30 Jun. 1986 (Illinois Univ.) 44 p CSCL 20L 

G3/76 


Prepared by 


S. L. Chuang 

Department of Electrical and Computer Engineering 
University of Illinois 
Urbana, IL 61801 

K. Hess, J. J. Coleman and J. P. Leburton 
Department of Electrical and Computer Engineering 
and Coordinated Science Laboratory 
University of Illinois 
Urbana, IL 61801 


3/VY3 


N87-1 2355 

Unclas 

44676 


SEMIANNUAL REPORT FOR 


Grant No. NAG 1-500 


SEMICONDUCTOR SUPERLATTICE PHOTODETECTORS 


January 1, 1986 to June 30, 1986 


Submitted to 


Dr. Ivan Clark 


National Aeronautics and Space Administration 
Langley Research Center, Mail Stop 473 
Hampton, VA 23665 


Prepared by 


S. L. Chuang 

Department of Electrical and Computer Engineering 
University of Illinois 
Urbana, IL 61801 

K. Hess, J. J. Coleman and J. P. Leburton 
Department of Electrical and Computer Engineering 
and Coordinated Science Laboratory 
University of Illinois 
Urbana, IL 61801 


2 


TABLE OF CONTENTS 


I. INTRODUCTION 

1. Period ......... 

2. Reporting Date . ...... 

3. Technical Personal ...... 

II. TECHNICAL PROGRESS 

III. FUTURE WORK 

IV. PUBLICATIONS 

V. TRAVEL 

APPENDIX A: Reprint of:” Impact Ionization Across the 
Conduction-Band-Edge Discontinuity of Quantum-Well 
Heterostructures" by S. L. Chuang and K. Hess, 

J. Appl. Phys., vol. 59, pp. 2885-2894, 1986. 

APPENDIX B: Copy of viewgraphs entitled: "Impact 
Ionization Across The Band-edge Discontinuity For 
A Superlattice Photomultiplier” by S. L. Chuang and 

K. Hess, presented at the 1986 Device Research 
Conference, Amherst, Massachusetts, June 21-23, 1986. 

APPENDIX C: Reprint of: "Lateral Shift Of An Optical 
Beam Due To Leaky Surface-plasmon Excitations" by 
S.L. Chuang , J. Opt. Soc. Am., vol. 3, pp. 593-599, 1986. 


page 

3 


4 

6 

7 

8 


CO CO CO 



3 


I. INTRODUCTION 

The research grant NAG 1-500 entitled "Semiconductor Superlattice 
Photodetectors" was awarded to the University of Illinois at Urbana-Champaign 
by National Aeronautics and Space Administration-Langley Research Center 
on June 27, 1984. The grant was continued on July 1, 1985. Dr. Ivan Clark 
is the Technical Officer, and Mr. John F. Royall is the Grants Officer. 

This report is the fourth semiannual report. 

1. Period: 

January 1, 1986 to June 30, 1986 

2. Reporting Date: 

July 1, 1986. 

3. Technical Personnel: 

S. L. Chuang Assistant Professor of Electrical and Computer 
Engineering 

J. J. Coleman Professor of Electrical and Computer Engineering 

K. Hess Professor of Electrical and Computer Engineering and 

Research Professor of Coordinated Science Laboratory 
J. P. Leburton Assistant Professor of Electrical and Computer 
Engineering and Research Assistant Professor of 
Coordinated Science Laboratory 


Two research assistants 


I 
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II. TECHNICAL PROGRESS 

During the past half-year period, our research effort centered on 
further investigating both the new superlattice photomultiplier with 
tunneling-assisted impact ionization and the photodetector based on the 
real space transfer mechanism. The following is a brief outline of the 
current status of these projects. Detailed analyses are included in 
Appendices A and B. New results of the tunneling-assisted impact ionization 
are presented in Appendix B and will be presented in a future publication. 


1 . Superlattice Photomultiplier 

We have, recently, included the effect of tunneling in the impact 
ionization across the conduction-band-edge discontinuity. The results 
have been presented at the 1986 Device Research Conference. A copy 
of the viewgraphs has been included in Appendix B. On pages B2-B4, we 
review the motivation for our work as well as some previous work. We then 
discuss our work on a new superlattice photomultiplier on pages B5-B7. 

The effect of tunneling is discussed on pages B8-B11. It is shown that, 
to properly take into account the tunneling process, a numerical approach 
using the Airy function, rather than the WKB method used commonly in the 
literature, should be followed because the conventional WKB method fails 
when the electron energy is near the barrier height. It is shown in page 
B-10 that, especially on the lower energy side, tunneling lowers the 
threshold energy and enhances the ionization rate. When the energy of the 
incident electron is very high, the final energies of the two electrons 
will be high, thus they tend to be above the barrier and can contribute 
to the current without tunneling. The multiplication factor versus the 
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applied voltage is shown on page B12. It shows a smoother increase than 
that of a general avalanche photodiode. Since our device involves only 
one type of carrier (electrons), it works as an ideal photomultiplier and 
does not have the problem of the rapid increase of the multiplication 
factor versus the applied voltage. The detailed theoretical formulation 
will be shown in a future publication. On pages B13-B14, the experimental 
results of Capasso et al . are shown to illustrate the plausibility of our 
predicted phenomenon; i.e., the impact ionization across the conduction- 
band-edge discontinuity. 

2 . Superlattice Photodetector Based on the Real Space Transfer Mechanism 

We are investigating the free carrier absorption process of a superlattice 
photodetector which makes use of the real space transfer mechanism. In 
particular, we have formulated and numerically computed the free carrier 
absorption coefficient for bulk GaAs in which we considered a second- 
order process involving both a photon and a phonon. Additionally, we 
carried out the computations for the free carrier absorption coefficient 
in a quantum well for the above cited second-order process involving a 
photon and a phonon. We are, currently, completing the numerical calculations 
for this case. Finally, we have begun to investigate the free carrier 
absorption coefficient in a quantum well in which we consider a first-order 
process involving only a photon. 



III. FUTURE WORK 


We shall apply our theory to the devices recently presented by 
Capasso et al . at AT&T Bell Laboratories. They have observed the predicted 
phenomenon of impact ionization in superlattice structures made of different 
materials, i.e., Al InAS/GalnAS , AlSb/GaSb, and InP/GalnAs. Though our 
theory has qualitatively predicted these phenomena, it will be interesting 
to compare the results quantitatively. The possible experimental approach 
will also be taken. 

In the area of the superlattice photodetector based on the real 
space transfer mechanism, we will continue to investigate the free carrier 
absorption process. The numerical computations of the free carrier absorption 
coefficient in a quantum well involving a second-order process of a 
photon and a phonon will be completed. The enhancement of the free carrier 
absorption coefficient for the quantum well over that of bulk GaAs will 
be examined. Additionally, we will continue our investigation of the free 
carrier absorption coefficient in a quantum well wherein a first-order 
process, employing only a photon, is involved. Consequently, we will be 
able to make a comparison between the free carrier absorption coefficients 
in a quantum well for the first- and second-order processes and determine 


which of the two is dominant. 
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IV. PUBLICATIONS 

The following manuscripts submitted for publication were supported 
either fully or partially by this Grant. The support by this Grant has 
been acknowledged in these manuscripts. 

1. S. L. Chuang and K. Hess, "Impact ionization across the conduction- 
band-edge discontinuity of quantum-well heterostructures," J. Appl . Phys . , 
vo 1 . 59, pp. 2885-2894, 1986. 

2. S. L. Chuang and K. Hess, "Impact ionization across the band-edge 
discontinuity for a superlattice photomultiplier," presented at the 1986 
Device Research Conference, Amherst, Massachusetts, June 21-23, 1986. 

3. S. L. Chuang, "Lateral shift of an optical beam due to leaky surface- 
plasmon excitations," J. Opt. Soc. Am., vol. 3, pp. 593-599, 1986. 

4. S. L. Chuang, "A coupled-mode formulation by reciprocity and a variational 
principle," IEEE J. Lightwave Technology, accepted for publication. 
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V. TRAVEL 

1. S. L. Chuang attended the Conference on Optical Fiber Communication 
and the topical meeting on Integrated and Guided Wave Optics, Atlanta, 
Georgia, February 23-28, 1986. 

2. S. L. Chuang presented a paper entitled: "Impact ionization across 
the band-edge discontinuity for a superlattice photomultiplier" at the 
1986 Device Research Conference, Amherst, Massachusetts, June 21-23, 
1986. He also visited AT&T Bell Laboratories in Murray Hill, N.J., June 
20, 1986. 
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Reprint of: "Impact Ionization Across the Conduction-Band-Edge 
Discontinuity of Quantum-Well Heterostructures" by S. L. Chuang and K. Hess, 
J. Appl . Phys., vol. 59, pp. 2885-2894, 1986. 


Impact ionization across the conduction-band-edge discontinuity of 
quantum-well heterostructures 

S. L Chuang 

Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, 
Illinois 61801 

K. Hess 

Department of Electrical and Computer Engineering and Coordinated Science Laboratory, University of 
Illinois at Urbana-Champaign, Urbana, Illinois 61801 

(Received 9 September 1985; accepted for publication 7 January 1986) 

Impact ionization across the band-edge discontinuity of quantum-well heterostructures is studied 
theoretically. We consider a heterolayer structure of alternating Al x Ga t _ x As and GaAs layers 
where the GaAs layers are heavily doped with donors. Thus a large number of electrons is 
confined to the quantum-well region. Incident electrons are heated up by applied electric fields 
and collide with the electrons confined in the well regions. Both the ionization rate as a function of 
the incident energy, and average ionization rates are computed. Device applications of such 
multiple quantum-well structures and the possibility of a complete analog to the conventional 
photomultiplier are discussed. 


I. INTRODUCTION 

Novel types of photodetectors using superlattice struc- 
tures have been proposed recently. 1-10 These include the en- 
hancement of impact ionization in a superlattice for ava- 
lanche photodiodes, 1-6 an infrared detector using free 
carrier absorption in a heavily doped quantum-well struc- 
ture, 7 a photodetector using the transient displacement cur- 
rent effect in a sawtooth superlattice, 8 strained-layer super- 
lattice photodetectors, 9 * 10 and doping (n-i-p-i) 
photodetectors. 11 * 12 The operational principles of all these 
new devices are very interesting and additional studies are 
needed to characterize these photodetectors completely. 

In this paper, we consider a heterolayer structure which 
consists of alternating Al* Ga { _ x As and GaAs layers where 
the GaAs layers are heavily doped with donors ( Fig. 1 ) . The 
donor are electrons confined mostly to the quantum wells. 
Incident electrons (which, for example, can be generated 
optically) are accelerated by the external field and collide 
with the electrons confined in the wells as shown in Fig. 1. 
This results in impact ionization of some of the electrons 
bound in the quantum wells. The wave functions for the 
incident electrons (three dimensional) and the confined 
electrons (quasi- two dimensional) are used to calculate the 
electron-electron interaction due to the Coulomb force with 
screening effects taken into account. The quantum mechani- 
cal transition rate for a system involving an initial state of 
one free electron above the well and one bound electron in 
the well and a final state of both electrons free is calculated 
using the Golden rule. The process is schematically dis- 
played in Fig. 2. The electron-electron interaction in the 
momentum space is illustrated in Fig. 3 and will be explained 
in the next section. The average ionization rate over the mo- 
mentum (energy) distribution of all the incident hot elec- 
trons is then calculated. Potential device applications are 
discussed in a final section. 

Our formulation for the impact ionization for the quan- 
tum-well structure is similar to the inverse Auger effect in 
bulk semiconductors 13-17 or in a quantum-well heterostruc- 


ture. 18-19 The new feature is that we are interested in impact 
ionization across the conduction-band-edge discontinuity 
while previous studies are considering transitions across the 
energy band gap. Thus, the threshold energy for impact ioni- 
zation is smaller (of the order of conduction-band-edge dis- 
continuity) than that of the interband process. The required 
bias voltage may be smaller since a smaller electric field than 
that of the conventional avalanche photodiode is needed. 

II. ANALYTICAL CONSIDERATIONS 

In this section, we present the detailed formulation of 
the impact ionization of electrons across the conduction- 
band-edge discontinuity in a heterolayer structure. First, the 
transition rate due to the electron-electron interaction is de- 
rived, then we discuss the energy distributions of the elec- 
trons confined in the quantum-well region and the incident 
hot electrons. Third, we present the impact ionization rate 
and its numerical evaluation. 

A. Transition rate 

Consider an incident free electron interacting with a 
bound electron in a quantum well; the total average transi- 
tion per unit volume P * is given by 19 * 20 



x[i-/(k;)][i-/(kj)], (l) 

where /is the Fermi distribution and P J is the quantum 
mechanical transition rate (per unit time) between the ini- 
tial state of electron 1 in state k, f electron 2 in state k 2 , and 
the final state of electron 1 in state k[ and electron 2 in state 
k 2 (Fig. 2) due to the electron-electron interaction: 

P?Z=2Z-\(l2\H,\V2mE l+ E 2 -E\ -E’ 2 ). (2) 

The square of the matrix element | < 1 2 \H S | IT > | ' 2 is given by 
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|<12|»,|1T)| ! 


-4 r«/r,)<(r,) ^-?l r '- r *'> 
J ^ 4ire |r, — r 2 | 


X ^ k j(r 1 )^ k ,(r 2 )rfr l ^r 2 , 


for the interaction by a screened Coulomb potential. Here q 
is the screening parameter which is discussed in Appendix 
A. We have neglected the difference in the dielectric con- 
stants of the two materials which is a good approximation 
for AlGaAs-GaAs. The factor of 4 in Eq. ( 3 ) takes account 
of effects of the exchange and preservation of the spin, 14 * 15 
considering the integrations over k( and k 2 , and 


FIG. 3. A geometrical configuration in k space to illustrate the electron- 
electron interaction integral. With both the energy and the momentum con- 
served, the four electron wave vectors are on the same spherical surface. The 
integrations over k| and k' 2 are restricted to the shaded portion with the z 
component between k t ' h and k lx 4- nrr/L « — k (n =* 1 ). 


<rir 2 |12>=^ k ,(r,)^ kj (r 2 ) 
_ exp(/k, • r,) 

W 


X exp(/ k 2 , • r 2 , )sin z 2 , 

“m 


<rir 2 |l'2*>=^ It , (r,)^. (r 2 ) 

_ exp(/k;-r,) exp(/k 2 • r 2 ) 

Jv 4V 


(4) 


(5) 


Initial State 



FIG. 2. Initial and final 
states for the electron-elec- 
tron interaction in a quan- 
tum-well structure. 



where k„ k [ , k 2f r lt r 2 are the general three-dimensional 
vectors and k 2f f r 2t are two-dimensional vectors in the quan- 
tum-well interface planes. The periodic parts n k (r) of the 
Bloch functions associated with Eqs. (4) and (5) have been 
ignored since the overlap integral is assumed to be one for 
states in the conduction band. The effective quantum-well 
width L zn is defined such that 



where E Qn is the energy level of the «th subband, which is 
calculated numerically taking into account the finite barrier 
height A E ct i.e. t the conduction-band-edge discontinuity 
(see Appendix B). For convenience, we further define 

k 2 * = k 2r ± z(nir/L zn ) = xk x + yk y ± z(nv/L zn ) , ( 7 ) 

r 2 = r 2l + zz 2 = xx 2 + yy 2 + zz 2 , (8) 

and replace the summations Z k by integrations. One then 
obtains 
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1 


( 9 ) 




and 


'■^Ttrl f <f> k -> J i£fc. nK > [S(l “’ k “ ,+S(k " k '- )1 ' 


( 10 ) 


Here we have used the fact that/(k; ) and/(kj )<1, also/Ck^ ) =/(k 2 + ) =/(k 2 - ) (which will be discussed later) and 


S(k„k 2 ±)s J 


dk| 8{E l + E 1 -E\ -Ei) 


(2ir) J 

^ tli 


212 


[(k^k tf + q 2 ) 


Ik2-k.+^ -k( 


(ID 


The lower limit in the integration for k{ is due to the fact that both k \ z and k ^ need to be larger than a certain threshold k ^ in 
order to get both electrons V and 2’ out of the quantum well. The condition k^k^ imposes another condition on k \ z 
because of the relation 


k 2r “ k\ z ± ~Z k lx ^ th * 

One finally obtains 

k’ th <k\ z <k lz ±nir/L zn -k' A . (12) 


One notes that the above condition for k \ z also holds for k ’ u obviously. 

An upper bound of the integral of Eq. ( 1 1 ) without the constraints on k \ z and k ^ has been obtained analytically for a 
parabolic band structure. 21 



dk; g(£,+£ 2 -£; -Ej) 
(2ir) 3 [|k,-k;| 2 + g 2 ] 2 


kj-k. + kjj* 


-k; 


m* 

8 v 2 #? 


Iki-kjSl 

|kj — k 2 ^ | 2 + ? 2 * 


(13) 


r 


The geometry pertinent to the evaluation of S(k,,k 2 ) 
and its upper bound S u (k„k 2 ) is shown in Fig. 3 for 

k 2 =k 2 t = k 2 , 4- z{mr/L zn ) (n = 1). 

Due to energy and momentum conservations, the four vec- 
tors k„ k 2 , k; , and k 2 lie on a spherical surface, assuming 
that the E — k relations are parabolic for all four vectors. 22 
Due to the above mentioned constraints on both k[ z and k^, 
the S(ki,k 2 ) integral will be over the unshaded spherical 
surface in Fig. 3 for which 

^ th 1 z ^ th 9 

while the integral for S u (k,,k 2 ) is over the whole spherical 
surface. The latter integration can be performed analytically 
by a coordinate transformation to the polar axis k l — k 2 , and 
results in Eq. (13). 


B. Electron distribution function f(k 7 ) in the quantum- 
well region 


The electron distribution function in the quantum-well 
region /( k^ ) or /( k 2 , ) is taken as 


AK) 

where 


[1 if0<* 2 ,<^‘,"> 
0 otherwise. 


(14) 


K iV = A2m*/#)E F - (nv/L z „ ) 2 (15) 

is the maximum of k 2 , for electrons in the nth subband. 
The electron concentration is given by 


n= f(E)g(E)dE 

JE<2> 

= \-^—k B T\a(l+e (E '- E °- )/k ‘ T ), (16) 

where E Qn is the /ith quantized energy level of electron 2 as 
given in Eq. ( 6 ) . Defining a characteristic concentration n c , 

n c =m*k B T/ir#L z . , (17) 

We have for E F — E Qn >k B T, 

— - X (18) 

n c occupied 
subbands n 

From n~N D and the above relation or Eq. ( 16), we obtain 
E Ff and thus K ( 2 *\ 

C. Electron distribution function /(k-i) for the incident 
hot electrons 

We assume that the incident hot electrons in the 
A l x Ga t _ x As region (region a) is Maxwellian in shape, i.e., 
the spherically symmetrical part is given by the expression 

/(ko) = 4« 0 ( rr tf/2m*k B T t ) 3,z e~ (19) 

with an electron temperature T e , where m* is the effec- 
tive mass of the electron in the conduction band in the 
Al x Gaj _ x As region and the index 0 characterizes the elec- 
trons in this region. Since the energy 
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E 0 = (#/2m»){k\ x +k 2 0y + k^) (20) 

is measured from the conduction-band edge in the 
Al* Ga, _ x As region, and 

V< ) = k 0x x + k 0y y + k (a z, (21) 

the energy E, and k, vector in the GaAs region (region b) 
are related to E 0 and ko by the equations 

( 22 ) 

=*o,» (23) 


tfn m*)kl + A E c = (#P/2 mt)k\, (24) 

where m% is the average effective mass of electrons in the T 
valley of the GaAs region. (We use the average effective 
mass m* to include the nonparabolicity effect. The effective 
mass m* used in the previous equations is replaced by m’ b in 
the computations. ( See Appendix C. ) Instead of integrating 
over d k„ we integrate over d ko, 




D. Average impact Ionization rate 

The average ionization rate ( 1/r) is defined as the total transition rate per incident electron, and therefore 

„ . C d\u .( it# ( . , 5 .\ 1 


< l/r> as Ptt/n 0 — f-^4f-^_y /J expf #- 

J (2 ir) 3 \2m*k B T,) V 2 m*k„T. 

= fef r af (- -***■-) ' 

1 2„ l 2 m‘k. T. ) V 2m**, T. ) r(S „ ) 


’ (k Ox + k 0y H" k o* ) 


where 




) — ! — 
J T(E 0 ) 


which is the ionization rate for an incident hot electron with 
the energy 

E 0 = E 0t = {ft L /2m*)kl t (27) 


= K = 0. (28) 

Since the integrand in Eq. (25) is sharply peaked near 
= koy —0, we have used the approximations 

— - — . , 0 = — 1 — , (29) 

r{E 0 ) t(E 0 ) r(£ te ) 

and carried out the dk 0x and dk 0y integrations in Eq. (25). 
By doing so, the incident electrons are essentially z directed. 

ko = zk Ql . (30) 

Thus, for the transmitted electron we have 

k { =zk lx (31) 

and 

(tf/lm^k^ + &E C = ( #/2m*)k \ x . (32) 

The average ionization rate ( 1/r) is obtained by integrating 
Eq. (25) using Eq. (26). An upper bound of <l/r>, (l/r) M , 
is obtained using the upper bound of \/r(E 0 ) . Dropping the 
subscript z and using E 0 instead of E 0x , one obtains for the 
integrand of Eq. ( 25 ) 

1 _ e 4 Air f^znY f dk 2t 

r(E 0 ) u e 2 ft\L) J 

X [S u (k|,k 2 „ ) + S u (k„k 2rt ) ] . (33) 

Next we first discuss the exact numerical integrations of 1/ 
r(£ 0 ) and (1/r) inEqs. (25) and (26). Then formulas for 
upper bounds, l/r(E 0 ) u , which lead to analytical expres- 
sions for ( l/r) u are presented. 


111. THE NUMERICAL INTEGRATIONS 

A. Integrations 1/r(E 0 ) and <1/r> from Eqs. (25) and (26) 

To obtain l/r(is 0 ), we need to evaluate the electron- 
electron interaction integrals S(k l9 k £ ) in Eq. ( 1 1 ). For the 
following discussion, we use k 2 =k 2 + . The integral 

S(k,,k 2 ~ ) is obtained similarly, replacing k 2 * by k 2 ~ , or mt/ 
L zn by — n Tr/L zn . For a parabolic band structure, and using 

k; = kj, +zk\ zy (34) 

K *kL+ifcL. (35) 

the 8(E X +E 2 — E\ —E\) term in the integrand of Eq. 
(11) can be rewritten as 

(«?/*)$[(*;, - *;?)<*'„ - *i ?)], (36) 

where the two roots for k \ z are 


k lz + mr/L z 


' k, z — nrr/L, 


k ,0 = 

K i* — 


k lz + nrr/L z 


’ k Xz — nir/L z 


-) + k !/ k 2, -*{?. 


■j 4- k lr • I 


Therefore, the k \ z integration can be performed analytically 
using the properties of the delta functions, and k \ z is re- 
stricted as given in Eq. (12). The other integrations can be 
simplifed as 
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J dk* j dk\, = J[* J ' kl, dk lt 

p K u r 2ir 

X Jo k\,dk\ t 2ir\ d<t> x 2, (39) 


where the two angular variables <f> lt for k 2t and <fr’ u for k|, 
have been transformed to only one angular variable <f> u 
which is the angle between kj, and k 2 ,, with a 2ir factor 
added as a result of this transformation because the inte- 
grand depends on <f> x2 only. The upper limit of k J, variable 
K [, is obtained from the condition that k J“ and k jf in Eqs. 
(37) and (38) should be real, i.e., 

: : I 


(— — nWL zn j + k j ^ cos ^ 12 _ k ;2>0. (40) 

This condition is easily included in the numerical code. 
If it is violated, k \ z and k f become inaginary and the inte- 
grand vanishes because the variable of integration k \ x is real. 
Explicit solutions of Eq. (40) can also be used, which give 
the upper limit for k J,, 




(41) 


Then l/r(£ 0 ) can be written as 


1 


e 4 mt 


t(E 0 ) e 2 # 


(L Y i r K *’ c 2w ( r K '* r 2 * 

J &2r dk 2t J d<f> X2 k lr dk [ f J d^ l2 F^(k\ tf k 2t9 ^i 2 ) 


+ JT k dk jf d<f>, 2 F ( k \„k ln <fi n ) j, 


(42) 


where 

F (k [ t9 k 2tt <f>\2 ) s 


l [(~ ^~~~~ ) + k U k 2. COS<fi t2 -k’J 


1/2 


([*;? 


1 


1 


'•+{k u -kS) 2 +f]* ’ [*{?+(* u -*i?) 2 + ^]V (43) 

for k \ z and between k ' h and k lz -4* (nir/L zn ) — £' h , and zero for the appropriate terms otherwise. F~(k [ n k 2n <fi l2 ) and 
K are obtained by replacing mr/L zn by — mr/L zn in the F * function as well as in K from Eqs. (41 ) and (43). 

The last integration J dk 0z in ( 1/r) must be performed from a minimum k ^ min which is obtained from Eq. ( 12) , 

k lz >2k' th ip nir/L zn , (44) 

and the energy conservation relation Eq. (32), 


k ± 

Oz,min 




(45) 


B. The upper bounds of the integrals 1/r(£ 0 ) u and <1/t>„ using Eqs. (25) and (33) 

The analytical result for the upper bound of 1/t(E 0 ) is 


1 


t(E 0 ) u 
G ± * 


e*m* b 

Yf? 


(¥)' 


i 


4 iPqL* 


(G + +G~y 


1 


1 



( k u Jf. nir/L zn \ 2 

yj + l 

\ q ) 


r(E 0 ) 
nir/L z 


t{E „) 


(46) 


J _ ^ k u ^ mr/L zn j 


— tan 


1 -^( A ’ 2 ') 2 1 ^ ku+ntr/Lu, J | tan" 1 j 


(47) 


Here the G + term is due to the term with , and the G 2 term due to that with Kn in Eq. (33), and 

r oc 

<l/r) u = 


r dk o. 

( IT# > 

1/2 

| ex Pl 

f #k i > 

1 1 

2tt 

U m*k a T.) 

l 2m*k B T, ) 

r(£ 0 ) 


r dkfa , 

f Y /2 / 

' # k l > 

1 1 


h^min 2ir ’ 

I 2 m*k B Tj P \ 

< 2 fn%k B T e ) 

1 r(E 0 ) 

u 


(48) 


which can be integrated numerically by Simpson’s rule. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 

The direct numerical integrations of the ionization rate 
l/r(£ 0 ) and the average ionization rate < 1/r) in Eq. (25) is 
straightforward. The upper bounds of the ionization rate 
and the exact numerically integrated results as functions of 
the energy of the incident hot electrons are shown in Figs. 
4(a) and 4(b) . One can see that the ionization rate increases 
rapidly when the energy E 0 is above the threshold value E, 




(49) 
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5000012nd Subbond 
3rd Subbond 



IONIZATION 
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ENERGY OF INCIDENT ELECTRONS (eV) 

FIG. 4. (a) Upper bound of the impact ionization rates l/r(£ 0 ) M as func- 
tions of the incident hot electron energy E 0 (m:E 0 , ) for each subband and 
the sum of all three subbands. The parameters used are L z — 200 A, 
L = 1200 A, N d = 5 X 10 18 cm -3 , x ~ 0.35, F=77K, T t = 2300 K (mj 
= 0.09m o , = 0.105™* = 0.262 eV). (b) Exact impact ionization 

rates from numerical integration of Eqs. (25) and (42). The same material 
parameters as those in (a) are used here. < 1 /r) u = 4.4 X 10 9 /s using Eq. 
(48"), and < 1/r) = 1.4xlOVs using Eqs. (25) and (42). 


for each subband n = 1,2,3- The overall minimum threshold 
energy for impact ionization will be due to the highest occu- 
pied subband of course. For a quantum-well size L z = 200 
A, A1 mole fraction x = 0.35 in the Al^Gai _ * As region, 
and N d — 5 X 10 18 cm" 3 , there are three occupied subbands 
among the six quantized levels in the well. The Fermi level is 
0.1092 eV measured from the conduction-band edge in the 
GaAs region. The threshold value for k ^ (the minimum 
value of k \ z or k ^ for electrons V and 2' to get out of the 
quantum well) is given by 

k ; h = yj2m£(AE c —8)/fi, ( 50) 

where <5 accounts for the tilting of the quantum well due to 
the large applied field (Fig. 1 ). Without space-charge effects 
we have 

S~eFL z , (51) 

where e is the magnitude of the electron charge. In all our 
computations we assume a field strength F = 10 kV/cm and 
a corresponding electron temperature of 2300 K at a lattice 
temperature of T =11 K. (No accurate data for T e in 
Al* Ga t _*As are available to our knowledge. The above 
value is a result of Monte Carlo simulations in GaAs. In Ref. 
23 there are some discussions on F-T e relations in 
Al JC Ga l As.) The necessary parameters used here are 
taken from Refs. 24 and 25. 

The average ionization rate < 1/r) involves the integra- 
tion over the energy E 0 and the distribution function. The 
upper bound of the impact ionization rate l/r(£ 0 ) tt as a 
function of the energy E 0 is shown in Fig. 4(a) for each 
subband and also for all subbands combined. These curves 
are much easier and faster to generate than the results of the 
exact numerical integration which are shown in Fig. 4(b) 
since the results in Eq. (46) are analytical expressions. They 
also can serve as a means of checking the numerical accuracy 
of the exact numerical integral in Eq. (42). This is because 
l/r(£ 0 ) tt fromEq. (46) is the upper bound of the numerical 
integral of l/r(£ 0 ) inEq. (42). l/r(£ 0 ) tt is the integration 
over the whole spherical surface in Fig. 3, for which all k\ 
and of the two electrons are allowed everywhere on the 
spherical surface. l/r(£ 0 )is the integration over only a por- 
tion of the spherical surface for which the condition 

^ th u ii flfr/ E zn — k t ' h 

is satisfied. If one relaxes this condition in Eq. (43), the 
results of integration (42) will be over the whole spherical 
surface, which is the upper bound. Thus a numerical test in 
Eq. (42) by simply removing the constraints on k [ z can be 
checked exactly with the analytical expression in Eq. (46). 
This process determines the number of discretized intervals 
for each integration variable in Eq. (42). One then adds the 
constraints on both k \ z and k ^ to obtain the exact numeri- 
cal integrals. 

A. The effect of multiple subbands 

As shown in Figs. 4(a) and 4(b), the ionization rates 
due to higher level subbands are clearly very important be- 
cause electrons in those subbands are in higher energy levels 
and thus easier ionized. Figure 4(a) shows the upper bounds 
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of the ionization rates for each subband and the total contri- 
bution is shown as the solid curve. The total rate in Fig. 4(a) 
has a discontinuity whenever either a new subband with a 
term corresponding to either + nir/L xn or — rtir/L m term 
comes in. The results for the exact numerical integrations for 
1 /t(E 0 ) are shown in Fig. 4(b) and the curves look smooth- 
er than those in Fig. 4(a). (The small fluctuations in the 
curves of Fig. 4(b) are due to numerical errors in the multi- 
ple integrations.) The electrons in the highest subband re- 
quire the lowest threshold energy for impact ionization and 
there are more incident hot electrons with enough energy to 
impact ionize these electrons. 

B. The effect of the quantum-well width L g 

In Figs. 5(a) and 5 (b), we decrease L z from 200 to 80 A 
keeping the doping concentration N D the same (5.0X 10 18 
cm" 3 ). The upper bounds of the ionized rates are shown in 
Fig. 5(a) and the exact numerical results are shown in Fig. 
5(b). There are three subbands for this structure and only 
the lowest one is occupied at this doping level. The ionization 
rate is smaller than the previous case in Fig. 4(b). Roughly 
speaking, it is because the number of electrons per unit area 
in the quantum- well region N D L Z is smaller when L x is de- 
creased and N d kept the same. Although the ground level of 
a thinner quantum well is higher than that of a thicker well in 
Fig. 4(b), the higher subbands in Fig. 4(b) also contribute 
to the total ionization rate and thus the ionization rate in the 
case in Fig. 4(b) is larger. 

C. The effect of the band-edge discontinuity 

In Figs. 6(a) and 6(b), we show the upper bounds of the 
ionization rates for all subbands and the results of the nu- 
merical integrations by decreasing the A1 mole fraction x 
from 0.35 to 0.25 with L z the same as that in Figs. 4(a) and 
4(b). We use the 60%-40% rule 26 ’ 27 for the conduction and 
valence-band-edge discontinuities in our calculation. If one 
uses a different rule with a number between 60% and 85% 
for the conduction-band-edge discontinuity, the ionization 
rate will decrease because the well is deeper, Once A E c is 
determined, the threshold energy for impact ionization can 
be obtained. Decreasing x will decrease A E CJ thus the 
threshold energy is lowered and the impact ionization rate is 
increased. This is clearly shown in Figs. 6(a) and 6(b) as 
compared to Figs. 4(a) and 4(b). 

O. The effect of doping and the total ionization rate 

If one increases the doping concentration in the quan- 
tum wells, one expects that the ionization rate will increase 
because there are more electrons available for impact ioniza- 
tion. This is shown in Fig. 7 where the dashed line is for the 
upper bound of the ionization rate and the solid line is from 
exact numerical integration. Increasing N D will increase the 
number of occupied subbands also since the Fermi level will 
be raised higher. For each N D , the computation takes about 
1500 cpu seconds as an average in the Cyber 175 Computer 
using about 1.6 X 10 6 number of discretized points (40 for 
variable k 2/ , 40 for k \ t , 20 for <f > , 2 , and 50 for k Qz ). When the 
Fermi level is raised, there will be more electrons in the high- 
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FIG. 5. (a) Upper bound of the impact ionization rate as a function of the 
incident hot electron energy E 0 for a different quantum-well width L x = 80 
A. (L «* 1200 A, x = 0.35, JV D =5Xl0“ cm" 3 , mj = 0.09m,,, m* 
— 0. 105 /Rq. A E e = 0.262 eV. ) There is only one occupied subband at this 
doping level, (b) Exact impact ionization rate as a function of the incident 
hot electron energy E 0 obtained from numerical integration of Eq. (25) 
using Eq. (42). The same parameters as those in (a) are used here. < 1/r) „ 
* 6.7 X lOVs, and <l/r> = 1.6 X 107s. 

er subbands also and they are easier to be ionized since the 
minimum threshold energy is decreased. The other effect is 
that a higher doping concentration may also increase the 
tunneling probability of the bound electrons. This may also 
enhance the ionization rate. 

V. POTENTIAL FOR DEVICE APPLICATIONS 

Avalanche photodioes which in their operational char- 
acteristic come close to the ideal of the conventional photo- 
multiplier have been proposed by Capasso 2 ^ 6 and have been 
termed solid-state photomultipliers. The principle of these 
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FIG. 6. (a) Upper bound of the impact ionization rates as functions of the 
incident electron energy E 0 for a smaller Barrier height {x — 0.25, 
A E c « 0.187 eV) than that in Fig. 4. We use m* = 0. 105m 0 , m? = 0.083 
m 0 for this value of x and A E c . The other parameters are L z — 200 A, 
L = 1200 A, N d = 5x 10 l8 cm -3 . (b) Impact ionization as functions of the 
incident electron energy E 0 . Same parameters are used as in (a). {\/r) u 
= 1.8X 10 ,o /s and <l/r> = 1.5X 10 9 /s. 


devices is based on a highly asymmetric ionization coeffi- 
cient for electrons (a) and holes ( 0). A conventional pho- 
tomultiplier does not, of course, involve any holes at all. The 
question arises therefore if the new technological possibili- 
ties which are offered by MOCVD and MBE crystal growth 
techniques can be used to construct the complete solid-state 
analogy of a photomultiplier. This question was one of the 
driving forces of our calculation in the previous sections. 
Insteady of photomultiplier electrodes one would use quan- 
tum wells and the secondary electrons would be emitted into 
the neighboring AlGaAs instead of into vacuum. The results 
for the ionization rate in Figs. 4-7 indeed show that appre- 
ciable gain can be achieved if the doping of the quantum well 



FIG. 7. Average ionization rate <l/r> (solid curve) and its upper bound 
< 1/r) „ (dashed curve) as functions of the doping concentration N D in the 
quantum well. We use L t — 200 A, L — 1200 A, x = 0.25, T =11 K, 
T t = 2300 K.The parameters m *, mj, and A£ c are the same as in Fig. 6(a). 


is — 10 19 cm“ 3 andthenumberoflayersis larger than ~ 100. 
From Fig. 7, the average ionization rate varies from 10 7 to 
10 10 1/s. Thus the electron ionization coefficients may vary 
from around 1 to 10 3 /cm or higher assuming the average 
electron velocity is 10 7 cm/s or lower. The ionization rate 
can also be increased by decreasing the band-edge discontin- 
uity. Primary electrons can be generated by the absorption of 
infrared radiation either by the electrons in the valence band 
(interband absorption as shown in Fig. 1) or by the quan- 
tum-well electrons which make a transition from bound to 
propagating states. The latter type of free carrier absorption 
is of first order (the quantum well provides the momentum 
and phonons are not necessary) and therefore significantly 
stronger than the free carrier absorption in bulk semicon- 
ductors. This fact has been noted by Chiu et a/. 28 who also 
show that multiquantum- well structures (similar to our 
structure in Fig. 1, except that they also dope the GaAlAs 
regions) can be used as infrared detectors. 7 Chiu et aL did 
not, however, investigate the possibility of gain by the ava- 
lanche process. But, the possibility of varying the band-edge 
discontinuity (and therefore the range of infrared detection 
almost at will in the range of energies 0<,ha) 50.4 eV) and 
the possibility of gain (analogous to the conventional photo- 
multiplier) are enormously attractive features. 

There are, however, also enormous technological prob- 
lems which need to be overcome to realize such a device. 
Since holes are not involved, the electrons which are ionized 
out of the well need to be replenished. This means that sepa- 
rate ohmic contacts have to be applied to a number of sets of 
quantum wells without changing the high resistivity of the 
AlGaAs layers in between. It may be possible to accomplish 
this complicated task with selective Schottky barrier con- 
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tacts which would be ohmic for the highly doped quantum 
wells. Another severe problem is the requirement of highly 
doped quantum wells neighboring to essentially depleted 
and highly resistive AlGaAs. Both problems are not basic in 
nature and do not seem out of reach considering the recent 
progress in microstructure technology. However, more 
studies (including also the effects of field ionization) are 
needed to confirm and show the device potential of such 
structures. 

Independently, Capasso has recently performed experi- 
ments and has obtained encouraging evidence for impact 
ionization of the kind discussed in this paper. These experi- 
ments may call for an extension of our theory to include 
effects such as the tunneling-supported impact ioniza- 
tion 29,30 due to the band bending as shown in Fig. 1. The 
bound electrons in the quantum well may acquire a small 
amount of energy by collision with an incident hot electron 
and subsequently tunnel into the conduction band in the 
AlGaAs region. This process reduces the ionization thresh- 
old energy in each subband. Thus the ionization rates, as, for 
example, shown in Fig. 4(b), may start from a lower thresh- 
old value for each subband and may also be enhanced be- 
cause of the tunnel-impact ionization process. More work 
needs to be done to predict the precise performance of such a 
structure. 

VI. CONCLUSIONS 

The impact ionization across the conduction-band-edge 
discontinuity of quantum-well heterostructures has been 
studied. A geometrical configuration to illustrate the contri- 
bution to the electron-electron interaction in the momen- 
tum space has been shown. The average ionization rate has 
also been calculated for various doping concentrations and 
geometries. We have included the contributions from the 
multisubbands and the nonparabolicity effects in our calcu- 
lations. The application of this impact ionization mechanism 
to a solid-state superlattice photomultiplier has also been 
discussed. 
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APPENDIX A: THE SCREENING PARAMETER g 

The dielectric screening parameter in the static limit is 
given by 31 

'■[tHI-SX' (a,) 

which is the Thomas-Fermi screening parameter, and it re- 
duces to l/( Debye length) for a Maxwellian-Bcitzmann 
distribution in the case of a nondegenerate electron gas. 

When the electron concentration is degenerate, and the 
temperature is low, we obtain the Thomas-Fermi screening 
parameter using a three-dimensional density of states: 

q= (m*e 2 /e&y l2 (3n/ir*) lf6 , (A2) 


where n is the electron concentration. 

If we consider a two-dimensional density of state 
g 2D CB)» we can cany out the integral in Eq. (Al) since 
#20 (£) xs a constant and obtain 

q=\m*f/e#) m {N/irL z ) m , (A3) 

where N is the number of occupied subbands. Our calcula- 
tions for a quantum well of size I, = 200 A and A£ c »0.3eV 
show that Eqs. (A2) and (A3) give almost the same result 
for a large n. 


APPENDIX B: THE QUANTIZED ENERGY LEVELSf*, IN 
A QUANTUM WELL INCLUDING THE EFFECTIVE-MASS 
DIFFERENCE 

The energy level Eq* is obtained by using the boundary 
conditions 


= fa, 

l _ l 
m* dz mf dz 


(Bl) 

(B2) 


at the two interfaces z = 0 and z = L t . The result is the well- 
known eigenequation: 

tan kL, = 2a/( 1 — a 2 ), (B3) 

where 

k=(2mfE) in /fi, (B4) 

a=[(m?/m*)(A£ c -£)/£] l/2 . (B5) 

The eigenequation (B3) includes both the even and odd so- 
lutions for the wave functions. The number of quantized lev- 
els in the quantum well is easily determined from the crite- 
rion. If 

(N-l)ir (2m?A£ c ) ,/2 L I Nir 

Z < < -T-t ("b) 

2 2 n 2 

the number of quantized levels = N. 


APPENDIX C: AVERAGE EFFECT MASS mg IN GaAs 
REGION (REGION b) WITH THE NONPARABOLICITY 
EFFECT INCLUDED 


The E-k relation including the nonparabolic effect is 
E(.l+aE)=tfk 2 /2m*, (Cl) 


where we use a = 0.67 (eV) -l , m* — 0.063 m 0 at 77 K, and 
m 0 is the electron mass in free space. 

Using the definition 


1 

m% (E) 


1 d 2 
# dk : 


(C2) 


and 





mJr (E)dE / AE C , 


(C3) 


where A E c is the conduction-band-edge discontinuity, we 
obtain analytically 


m? = m*(l + aA£ c .)(l +2a&E c + 2a 2 A£ c 2 ). (C4) 


The above procedure accounts for the nonparabolic effect up 
to the energy A E c . We then use this average value mf in the 
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parabolic band structure for GaAs in our calculation for the 
ionization rates. 
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Copy of viewgraphs entitled: "Impact Ionization Across the 
Discontinuity For A Superlattice Photomultiplier" by S. L. 
K. Hess, presented at the 1986 Device Research Conference, 
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Conclusions: 


/g/s" 


A new superlattice photodetector has 
been proposed and studied theoretically. 

A single carrier type (electron) is involved 
and it has the potential of minimizing the 
excess noise and achieving an appreciable 
gain. 
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A leaky-wave theory is applied to the study of the lateral shift of a bounded optical beam reflected from a 
priam-air-metal (and prum-metal-eir) layered configuration with surface-plasmon excitations. The locations of 
the pole and the zero of the reflection coefficient for the leaky surface-plasmon mode are calculated as one varies the 
thickness of the center layer. The excited surface-plasmon mode is leaky because of the presence of the prism. 
Both the real part and the imaginary part of the pole are affected by the energy leakage. The detailed structure of 
the reflected intensity and the forward and backward shifts of the reflected beam are illustrated. 


1. INTRODUCTION 

Since the optical properties of a metallic film or a semicon- 
ductor sample can be characterized by a permittivity func- 
tion «(w), the theories and techniques to determine <(«) have 
been important to our understanding of the optical proper- 
ties of solids. One useful optical technique is to use the 
surface-plasmon wave theory in the reflectance measure- 
ment. 1 The excitation of surface-plasmon waves is achieved 
either by using the attenuated-total-reflection method or by 
using a grating structure with a proper choice of the phase- 
matching condition. If one measures the reflected power as 
a function of the incidence angle or the frequency, anoma- 
lous absorption will occur when the surface-plasmon mode is 
excited. Thus the surface-plasmon theory provides a sensi- 
tive tool to determine the permittivity function of the solids. 
For metals, the surface-plasmon mode is excited in the visi- 
ble frequency region, whereas for degenerate semiconduc- 
tors it is excited in the infrared region. The study of sur- 
face-plasmon excitations may give detailed information 
about the optical properties of the metals and semiconduc- 
tors. 

However, so far, almost ail the studies on the excitation of 
the surface-plasmon theory assume that the incident wave is 
an infinite plane wave. In electromagnetics or in acoustics, 
it is known that lateral beam displacement may occur when 
the incident wave is a bounded beam and proper phase- 
matching conditions are satisfied. 2 For example, if the inci- 
dent angle is at the critical angle between two media, the 
well-known Goos-Hanchen shift occurs. 3- * In acoustics, if 
the incident acoustic beam is from a fluid medium to a solid 
surface and the angle of incidence is near the Rayleigh criti- 
cal angle, the excitation of the Rayleigh wave in the solid, 
which becomes a leaky wave because of the existence of the 
fluid, causes a beam shift in the reflected wave. 7 The result 
of the reflected field is different from that predicted by the 
geometric acoustics. The intensity and phase variations of 
the shifted beam give information about the acoustic prop- 
erties of the solids. Thus this idea can be applied in the 
nondestructive testing of the mechanical properties of sol- 
ids. A leaky-wave theory has been successful in explaining 
various observed phenomena in acoustics 7 and for multilay- 
ered or grating structures. 8-10 


The properties of surface polaritons in layered structures 
have, been discussed in many references. 1 - 11 * 13 The attenu- 
ated-total-reflection method 1 is used to determine the opti- 
cal constants of thin metallic Aims 14 or a CsBr sample in the 
infrared region. 15 Grating structures have also been used to 
study the dispersions of surface plasmons in semiconduc- 
tors 16 - 17 and in metals. 16 Both prisms and gratings have 
been used to couple the electromagnetic wave beams into the 
waveguide modes. 19-21 Coupling of the finite beam with the 
surface-plasmon mode has also been studied, including the 
coupling efficiency 22-24 and the beam-shift phenomenon. 25 - 20 
The study on the lateral beam displacement in Ref. 26 close- 
ly follows Ref. 6, which assumes a Taylor-series expansion 
near the incident angle and is applied to the case of the 
critical angle of incidence, i.e. f the Goos-Hanchen shift. 
That expansion is not valid if there is a pole, which is due to 
the leaky surface-plasmon mode, as will be explained later, 
near the horizontal wave number where the reflection coeffi- 
cient is expanded into the Taylor series. This has been 
pointed out in Ref. 7, and a better approximation, which 
takes into account the singularity of the reflection coeffi- 
cient, has been proposed. The pole of the reflection coeffi- 
cient accounts for the leaky surface mode excited by the 
incident beam. The excitation of the surface electromag- 
netic waves in an Otto configuration (i.e., prism-air-metal) 
by using the leaky- wave theory has been studied 25 both 
theoretically and experimentally. In this paper, we apply 
the leaky-wave theory to both the Otto configuration and 
the Kretschmann configurations (i.e., prism- metal-air). 
The trajectories of the poles and zeros as the middle-layer 
thickness is varied are studied in detail. Other singularities, 
e.g., the branch points, are also shown. When the branch- 
point contribution is important in certain cases of Kretsch- 
mann configuration for which the leaky-wave theory for the 
reflection coefficient is no longer valid, a direct numerical 
approach is adopted, and the results are also shown. 

2. FORMULATION 

The theoretical formulation closely follows Ref. 10, except 
that we deal with the excitation of a leaky surface-plasmon 
mode that requires the incident wave to be TM polarized 
and either medium 2 or medium 3 has a negative permittivi- 


0740-3232/86/050593-07$02-00 •<* 1986 Optical Society of America 



594 J. Opt Soc. Am. A/VoL 3, No. 5/May 1986 


9RI6INAW PAGE is 
POOR OuALltY 


S. L. Chuang 



Fi«. 1- The reflection of a bounded beam from a priam-air-metal 
or pmm-metal-air configuration. 


ty in the real part that is smaller than -<o- The geometrical 
configuration of the problem is shown in Fig. 1, where both 
prism-air-metal and prism-metal-air configurations are 
considered. A metal such as silver has a permittivity ( - 5. 19 
+ i0.28) eo at the wavelength 4358 A. 14 

Assume that a bounded beam is incident from the prism 
region with the magnetic field given by 

H, - yH t (x, 0) « y exp +i *i i x J ( l > 

on the z » 0 plane, where W x » W/ coe 0„ W is the beam- 
width, &i is the angle o f inci dence for the center of the beam, 
k xi - ki sin 0„ kx » and is the permittivity of the 

prism. 

Using the plane-wave spectrum of the incident field at z * 

0 

*,(*, <»-{'£ h(k x )np(ik x x), (2) 

where 

h(k x ) - &W, e*p[-(*, - Kf O) 

the reflected magnetic field is given by 
H r « yft r (x, z) - y j R(At I )/»(fc I )exp{i* r i + ik ir z], 

(4) 

where k l2 - («Vo«i “ k x 2 ) in , with Im (k u ) i 0, and the path 
of integration is on the real axis of the complex k x plane. 
R(k x ) is the plane- wave reflection coefficient from the lay- 
ered medium. From the expression in Eq. (3), one sees that 
the major contribution to the integral in Eq. (4) is from the 
visible region with k x satisfying 

(5) 

As will be explained in Section 3, if 0* is near the angle 
where the surface plasmon is excited, the reflection coeffi- 
cient can be expressed as the product of a factor containing a 
zero and a pole in the form 


MJaRtr—r' (6) 

R x — R p 

where Rq is independent of k x as an approximation, provided 
that there are no other singularities such as the branch 
points in the visible region. If there are branch points near 
the region in expression (5), the approximation in expression 
(6) is no longer valid, and a direct numerical approach using 
Eq. (4) should be taken. This will be the case for the 
Kretschmann attenuated-total-reflection configuration, as 
will be shown later. The zero and the pole of R(k s ) in the 
complex plane can be obtained numercially, since an exact 
analytical expression of R(k x ) is known. Ro can be obtained 
accurately from R(A X i)(*xi - k p )/(k xi - k Let us denote 
k p and k n as 

fcp * + ict p , (7) 

k n m + »'«»• W 


The reflected field at z * 0 can be evaluated from Eq. (4) and 
expression (6) to be 10 

H r (x, 0) - R 0 + i & (k p - k n )W x exp<T 2 )erfc<T) 

X CXP [ - (^) + ****]’ (9> 

where 


7 



+ -0p> IT 


(10) 


and erfc( 7 > is the complementary error function of complex 
' argument 7- 


3. THE PLANE-WAVE REFLECTION 
COEFFICIENTS FROM A PRISM-AIR-SILVER 
AND A PRISM-SILVER-AIR CONFIGURATION 

The plane-wave reflection coefficient R(k x ) from a layered 
medium shown in Fig. 1 can be written as 27 

R \2 4 “ ^12^23 exp(i 2 Afl,d) 

— — - — ~ (id 

R 12 ^23 4 " Ru ^23 exp(i2A 2l d) 

where 

kj- k iT 

R,j* a — ± — . (12) 

‘i «> 

*„=.(« Vi -*x 2 ) 1/2 - (13) 

One notes that there are only two branch points in Eq. (11) 
occurring at k u * 0 and k & =■ 0. This has been noted in Ref. 
10 and also proved by invoking the uniqueness principle in 
solving the boundary- value problem for stratified media. 28 
Another algebraic method is simply to look at Eq. (11) and 
change k<u to one obtains the same expression for the 
reflection coefficient. Let us define the normalized quanti- 
ty 

*~k x /k, (14) 

where k * is the wave number in free space. The two 

branch points occur at 
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Fig. 2. (a) The branch cuta and the location! of the branch point!, the pole kp/k f and the zero hjk in the complex « plane for a prism- 
air-silver configuration, (b) The trajectories of the poles and zeros on the xjriane when the air-gap thickness is varied. 
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Fig. 3. (a) The branch cuts and the locations of the branch points, the pole k p /k, and the zero kjk in the complex * plane for a prism-silver-air 
configuration, (b) The trajectories of the poles and the zeros on the complex * plane when the silver thickness is varied. 


* V * 1^*0 ( 15 ) 

and 

^3 s V *o (15) 

in the complex * plane. The branch cuts are chosen to be 27 
Re[A u ] -0 (17) 

and 

ImlAJ - 0, (18) 

which are shown in Figs. 2(a) and 3(a). 

Case 1. Prism-Air-Metal (Otto Attenuated Total 
Reflection) Configuration 

The branch cuts, as defined by Eqs. ( 17) and (18), are shown 
in Fig. 2(a). Once these branch cuts are defined, the pole 


and the zero of the reflection coefficient for k x in Eq. (11) can 
be obtained and written as kp/k and kjk in the « plane after 
being normalized to the free -space wave number k. 

The pole in the reflection coefficient k p corresponds to a 
leaky surface-plasmon mode of this configuration. If the 
prism is not present, the surface-plasmon mode has a disper- 
sion relation 



If both <2 And <3 are reai and « 3 < -< 2 » the above solution in 
Eq. (19) is real and corresponds to a nondissipative surface- 
plasmon mode. Adding the prism will make this pole com- 
plex because it corresponds to a leaky mode, i.e., the cou- 
pling of the prism and the air-silver surface-plasmon mode 
would make this mode radiate into the prism region. A 
complex pole comes from the fact that it is a leaky mode but 
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not from a power dissipation in the bulk region because <i 
and <3 are real. 

If <3 ie a complex quantity, there is a dissipation in the 
silver region, and the air-silver surface-plasmon mode with- 
out the prism has (kjk)&^ 9 m 1.1125 + <0.0071 from Eq. 
(19). With the prism, the leaky surface-plasmon mode that 
has the wave number kjk close to contains both 

the leakage and the dissipation effects. As d/X— «•, we have 
exp(i2*2rd) 0 [one chooses Im(At2*) > 0 for convenience] 
and 


*(*,)■ 


^ 12 ^ 23 * 


Mr 


' 12 * 


( 20 ) 


Thus both the pole and the zero of R(k x ) approach the root of 
#23* “ 0 that gives exactly Eq. (19) as it should, since, as d » 
X, the reflection coefficient R(k x ) should be independent of 
the property of the air-silver interface, and the pole and the 
zero must cancel each other. On the other hand, if d -* 0, 
the reflection coefficient can be rewritten as 


#(* x ) 


<1 *3 

-k u + ±k„ 

‘l <3 




( 21 ) 


which has a pole at flu* ■ 0, Le., 



and we obtain kp/k* 2.0604 + <0.0454. Thus the pole in Fig. 
2(b) will approach this value as d/X — 0. The zero of the 
reflection coefficient does not exist on the complex plane 
shown in Fig. 2(b). It actually moves to the bottom Rie- 
mann sheet defined by Re[&u] < 0 (note that we always 
define Im(k&) £ 0 in our discussions here] as it crosses the 
branch cut defined by Re{k u ) 9 0 shown in Fig. 2(a). This 
zero will then approach the value 2.0604 + <0.00454 in the 
bottom Riemann sheet with Reiku) < 0. 


Case 2. Prism-Metal-Air (Kretschmann Attenuated 
Total Reflection) Configuration 

The branch cuts are shown in Fig. 3(a). On the complex « 
plane, the branch cut Imffcs,] - 0 extends from the positive «/ 
axis to the positive axis with 0 < £ n$ » 1. The other 

branch cut Re(k u ) * 0 extends from ** * n x to on the 
real axis since <i is reaL The locations of the pole and the 
zero move as one varies the thickness of the metal d. As d/X 
— ® [Fig. 3(b)], both the pole and the zero move to the value 
defined by Eq. (19), i.e., ( kjk ) Ag _tir * 1.1125 + i0.0071, as 
can be seen from Eq. (20) with the common factor #23* * 0 
giving the pole and the zero in the range k t of our interest 
[expression (5)]. As d/X — 0, #(A X ) approaches #i3~/#i3* as 
in Eq. (21), and we have 



j - 0.8370 - <5, (24) 

where 5 is an infinitesimal positive number. In this case, the 
pole and the zero do not cancel each other because one is 


slightly above and the other slightly below the branch cut 
Im(*3t) ■ 0. 

Another difference between the Kretschmann configura- 
tion and the Otto configuration is the branch-point contri- 
bution. Aa shown in Fig. 3(b) for d/X between 0.02 and 0.03, 
the pole has a real part cloee to the branch point /13 * 1.0. 
Thus, if the angle of incidence is such that is close to the 
branch point, the approximate formula in expression (6) for 
the reflection coefficient will not be valid since the branch- 
point singularity will also affect the integral for the reflected 
field. In this case a direct numerical integration using Eq. 
(4) is necessary, and the results are also given in Section 4. 
Both a lateral-wave contribution (due to the branch point) 
and the surface-plasmon-mode contribution (due to the 
pole) are important in determining the reflected beam. 

4. NUMERICAL RESULTS AND DISCUSSIONS 

In Figs. 4-6, we show the reflected-beam intensities |H P (x, 
0)| 2 (solid curves) and the incident-beam intensities | Hi(x, 
0)|* (dashed curves) as functions of the position x for the 

prism l 

Oir d 

silver J" 





Fig. 4. Optical beam shift for the prism-air-silver configuration 
with d/X- 0.4 for four beam widths. The dashed curves represent 
the incident-beam intensities | H t {x, 0)| a , and the solid curves are for 
the reflected-beam intensities \H r {x, 0)j 2 . 
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Fig. 7. Optical beam shift for the prism-siiver-air configuration 
with WJ\ * 20 and four different silver-film thicknesses, dl\ - 
0.01, 0.08, 0.15, and 0.2. 
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Fig. 8. Optical beam shift for the priam-ailver-air configuration 
with dl\ * 0.025. Aa can be seen from Fig. 3(b), the location of the 
pole ia right above the branch point at « » 1.0. Thus a direct 
numerical integration of Eq. (4) ia required to obtain the reflected 
field. Dashed curves, incident- beam intensities; dotted curves, re- 
fiected-beam intensities using Eq. (9); solid curves, reflected- beam 
intensities using the direct numerical integration. 


tion, as demonstrated both theoretically and experimentally 
at microwave frequency in Ref. 9. 

In Figs. 7 and 8, the priam-silver-air configuration is con- 
sidered. D7A for both forward and backward beam shifts 
are calculated in the table of Fig. 3(b). The detailed intensi- 
ty variations are illustrated in Fig. 7 for different metal 
thicknesses, d/\ » 0.01, 0.08, 0.15, and 0.2. When d/X » 
0.025, the branch-point contribution from « - n 3 » 1.0 ia 
important As mentioned in Section 3, the approximation 
for the reflection coefficient in expression (6) in the leaky- 
wave theory 7 ' 10 is not valid. A direct numerical integration 
using Simpson’s rule for the reflected field is performed by 
using the exact reflection formulas in Eqs. (4) and (11). The 
beam shift still occurs as expected. This case is different 
from the Gooe-Hanchen shift in the contributions from both 
the lateral wave (branch-point contribution) and the sur- 
face-plasmon mode (pole contribution) in our case. The 
Gooe-Hanchen shift comes only from the branch-point con- 
tribution. If one still uses the simple pole approximation, 
the results (shown as the dotted lines) are quite different 
from the exact numerical solution (shown aa the solid lines) 
in Fig. 8. For the case in Fig. 8, the simple leaky-wave 
theory tends to overestimate the reflection coefficient and 
underestimate the beam shift, if the branch-point contribu- 
tion is ignored. 

At other wavelengths, e.g., X - 6328 A, « (-16.32 + 
t0.54) cq and k 9p /k at (1.0321 + iO.OOll). Thus the leaky 
surface- plasmon pole will be very close to the branch point 
at n 3 ■ 1.0 for the Kreschmann configuration. Again the 
leaky-wave approximation in expression (6) will not be ap- 
propriate, and a direct numerical integration, as has been 
done for Fig. 8, is necessary. 

5. CONCLUSIONS 

The leaky-wave theory has been applied to study both the 
Otto and the Kreschmann attenuated-total-internal-reflec- 
tion configurations, when the leaky surface-plasmon modes 
are excited. It is demonstrated theoretically that the beam 
shift and intensity variation due to surface-plasmon excita- 
tion exist, as predicted from the leaky-wave theory. This 
beam-shift phenomenon may affect the reflectance mea- 
surement if the beam width is narrow (for example, much 
less than 100 wavelengths). A study of the beam-shift phe- 
nomenon may also provide information about the dielectric 
functions of the metal or the semiconductor in the visible or 
the infrared region, respectively. In addition to the prism 
configurations, grating structures may also be employed for 
similar studies. 

ACKNOWLEDGMENTS 

Stimulating and helpful discussions with G. A. Des champs 
are gratefully appreciated. This research was partially sup- 
ported by grants NAG 1-500 from the National Aeronautics 
and Space Administration and ECS 83-11345 from the Na- 
tional Science Foundation. 


REFERENCES 

1. A. Otto, in Optical Properties of Solids: New Development , B. 
O. Seraphin, ed. (North- Holland, Amsterdam, 1975). 

2. R. Tamir, Integrated Optics , 2nd ed. (Springer-Verlag, New 
York, 1979), Chap. 3, pp. 83-137. 


4 


I 

i 


S. L. Chuang 

3. Von F. Goos and H. Hanchen, “Ein neuer und fundamentalcr 
Versuch zur Totalreflexion,” Ann. Physik 6, 333-346 (1947). 

4. H. K. V. Lotsch. “Reflection and refraction of a beam of light at 
a plane interface/* J. Opt. Soc. Am. 56, 551-561 (1968). 

5. T. Tamir and A. A. Oliner, “Role of the lateral wave in total 
reflection of light/’ J. Opt. Soc. Am. 59, 942 ( 1969). 

6. B. R. Horowitz and T. Tamir, “Lateral displacement of a light 
beam at a dielectric interface/’ J. Opt. Soc. Am. 61, 586-594 
(1971). 

7. H. L. Bertoni and T. Tamir, “Unified theory of Rayleigh -angle 
phenomena for acoustic beams at liquid-solid interfaces/* Appl. 
Phys. 2, 157-172(1973). 

8. T. Tamir and H. L. Bertoni, “Lateral displacement of optical 
beams at multilayered and periodic structures,” J. Opt. Soc. 
Am. 61,1397-1413(1971). 

9. A. Amittay, P. D. Einziger, and T. Tamir. “Experimental obser- 
vation of anomalous electromagnetic absorption in thin- layered 
media/’ Appl. Phys. Lett. 38, 754-756 ( 1981). 

10. V. Shah and T. Tamir, “Absorption and lateral shift of beams 
incident upon lossy multilayered media,” J. Opt. Soc. Am. 73, 
37-44 (1983). 

11. P. Halevi, “Polariton modes at the interface between two con- 
ducting or dielectric media,” Surface Sci. 76, 64-89 (1978). 

12. D. L. Mills and A. A. Maradudin, “Properties of surface polari- 
tons in layered structures,” Phys. Rev. Lett. 31, 372-375 ( 1973). 

13. A. D. Boardman, ed.. Electromagnetic Surface Modes (Wiley, 
New York, 1982). 

14. W. P. Chen and J. M. Chen, “Use of surface plasma waves for 
determination of the thickness and optical constants of thin 
metallic films.” J. Opt. Soc. Am. 71, 189-191 ( 1981). 

15. M. Tacke, “Far infrared laser spectroscopy of surface (Fano) 
waves for the determination of the optical constants of solid,” 
Sov. J. Quantum Electron. 12, 662-664 (1982). 

16. M. MarschalL B. Fischer, and H. J. Queisser, “Dispersion of 
surface plasmons in InSb,” Phys. Rev. Lett. 27, 95-97 (1971). 


Vol. 3, No. 5/May 1986/J. Opt. Soc. Am. A 599 

17. P. M. van den Berg and J. C. M. Bor burgh, “Dispersion of 
surface plasmons in InSb-gratings,” Appl. Phys. 3, 55-60 
(1974). 

18. R. H. Ritchie, E. T. Arakawa, J. J. Cowan, and R. N. Hamm, 
“Surface -plasmon resonance effect in grating diffraction,” 
Phys. Rev. Lett. 21, 1530-1533 (1968). 

19. P. K. Tien and R. Ulrich, “Theory of prism-film coupler and 
thin-film light guides,” J. Opt. Soc. Am. 60, 1325-1337 (1970). 

20. R. Ulrich, “Theory of the prism-film coupler by plane-wave 
analysis,” J. Opt. Soc. Am. 60* 1337-1350 (1970). 

21. M. L. Dakss, L. Kuhn, P. F. Heidrich, and B. A. Scott, “Grating 
coupler for efficient excitation of optical guided waves in thin 
films,” Appl. Phys. Lett. 16, 523-525 (1970). 

22. D. L. Begley, D. A. Bryan, R. W. Alexander, and R. J. Bell, 
“Optimization of prism coupling efficiency for surface electro- 
magnetic wave excitation in the infrared,” Appl. Opt. 16, 
1549-1552(1977). 

23. R. Miller, D. L- Begley, and R. J. Bell, “Surface electromagnetic 
wave coupling efficiency dependence on prism gap height,” 
Appl. Opt. 16, 3077-3079 (1977). 

24. R. T. Deck, D. Sarid, G. A. Olson, and J. M. Elson, “Coupling 
between finite electromagnetic beam and long-range surface- 
plasmon mode,” Appl. Opt. 22, 3397-3405 (1983). 

25. W. P. Chen. G. Ritchie, and E. Burstein, “Excitation of surface 
electromagnetic waves in attenuated total- reflection prism con- 
figurations,” Phys. Rev. Lett. 37, 993-997 (1976). 

26. P. Mazur and B. Djafari-Rouhani, “Effect of surface polaritons 
on the lateral displacement of a light beam at a dielectric inter- 
face,” Phys. Rev. B 30, 6759-6762 (1984). 

27. J. A. Kong, Theory of Electromagnetic Waves (Wiley, New 
York, 1975), especially p. 119 for reflection coefficient and p. 
213 for branch cuts. 

28. W. C. Chew, “The singularities of a Fourier- type integral in a 
multi cylindrical layer problem,” IEEE Trans. Antennas Prop- 
ag. AP-31, 653-655 (1983). 


